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Abstract 

We find a new stable rotator solution in DO-brane matrix mechanics. The 
solution is interpreted as a D2/D0 brane bound state, constructed as a trans- 
versely rotating ellipsoidal membrane with N DO-branes pinned over it. From 
the membrane point of view the attractive force of tension is exactly can- 
celled by the repulsive centrifugal force. Dynamical properties of the system 
are investigated, in particular the emission of spherical waves in the different 
massless sectors of supergravity. We compute the RR 1-form quadrupole, the 
RR 3-form dipole and the gravitational quadrupole radiation. Also, we show 
that our non-singular classical solution is stable against small perturbations 
in the initial conditions. Furthermore, we comment on a possible fundamental 
particle interpretation of our system. 
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1 Introduction 

We find a new stable D2-D0 brane bound state of the non-commutative fuzzy 
sphere. This is realized as an ellipsoid with three independent angular mo- 
menta corresponding to rotation in three different planes. The idea is to cancel 
the attractive force due to tension with the repulsive centrifugal force. As a 
result, the three angular velocities are uniquely determined through the princi- 
pal radii of the ellipsoid. We compare this to the spherical membrane solution 
|l[ where the sphere oscillates in a manner that inverts the orientation of the 
sphere at the point where the geometry is singular. In contrast, our solution 
is non-singular at all times. 

Before presenting our new solution in Section |3] we give in Section ^J] the 
necessary background on the matrix mechanics of DO-branes, and in particular 
the new approaches [fj, |3|] to the coupling of DO-branes to higher form RR 
fields. It is sufficient to consider the 3-form RR field which corresponds to D2- 
brane charge. In Section |2.2j , we review the spherical membrane solution, the 
essential elements of which we use in our construction of the rotating ellipsoidal 
membrane in Section [3[ 

In Section [| we analyze the classical stability of our solution against small 
perturbations in the initial conditions. It is known that the diagonal oscillating 
sphere solution is unstable [|], [J, resulting in the system spending most of the 
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time in long, string-like configurations. In contrast, our rotating solution is 
stable against small perturbations of initial conditions, as evidenced by the 
absence of non-unitary Poincare-Lyapunov exponents in the S-matrix. We use 
the full machinery of the Hamiltonian theory of dynamical systems to isolate 
the zero modes of the variational equations of motions. All the remaining 
eigenvalues were measured numerically and found to be unitary. 

The system exhibits interesting dynamical properties due to the interaction 
with various supergravity fields. This leads to semi-classical energy loss, and 
we compute in Section [5] the functional form and total power emitted by the 
system in three different sectors. 



Section |5.1| is devoted to discussion of radiation carried by the RR 1-form 
field C^. The following phenomenon appears: even though DO-branes nat- 
urally couple to this field and create a static potential, the dipole moment 
vanishes because only one kind of charge is present. Therefore we look at 
quadrupole order in search of radiation, and find that for the case of a 

sphere of radius R rotating at the frequency lu, the power is proportional to 
12 



2 8 2k 2 N 2 T? 



u 12 R 4 . (V 
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In Section |5.2| we use the machinery developed in Section |2.1| , and consider 
the 3-form RR field radiation which corresponds to D2-brane charge. A spher- 
ical membrane would not carry a net charge due to cancellation of opposite 
oriented pieces, but a dipole moment is non- vanishing in accordance with naive 
expectation. This rotating dipole produces radiation which in most respects 
behaves like ordinary dipole radiation, with only difference being that the 
gauge field has three indices instead of just one, in the case of a vector gauge 
field. This broadly fits with the understanding of RR gauge fields as corre- 
sponding to photon-like "spin-one" particle, even though there is no invariant 
meaning to the notion of spin in 10 dimensions. The power is proportional 
to u> 10 , but owing to the specific relation between R and u for our system 
turns out to differ only in the numerical coefficient 

o7 n 2 N 2 T 2 

P -2W^^ aR ■ (2) 



For gravitational radiation we find in Section p.3| radiation at the lowest 
possible order, namely quadrupole. The final result is 

2 6 2k 2 N 2 T 2 12 4 

°V 2 i? 4 . (3) 
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We note that the remaining massless field, the dilaton, at lowest order 
couples to F 2 and this quantity turns out to be time- independent. However, 
at higher order the dilaton may couple to time-dependent combinations of F. 
The calculation would require knowledge of higher order terms in the non- 
abelian DBI action, and we leave this to future work. 
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Following (El, we discuss the fundamental particle interpretation of our 
system in Section It is possible that there exists a completely stable ground 
state, protected by certain conservation laws from further decay This particle 
would appear to the external observer to have N units of DO-brane charge, and 
(presumably quantized) internal angular momentum. 



2 Preliminaries 



2.1 The dynamics of N DO-branes 

The system of N interacting DO-branes in type IIA string theory has been 
studied extensively in recent years. One of the reasons for this is the conjecture 
0, H that it should describe the discrete light-cone quantized (DLCQ) M- 
theory. 

Before presenting the action for N DO-branes we first shortly review the 
dynamics of the background fields. The background fields for the N DO-branes 
are governed by the type IIA low energy effective action 



Sua = —2 I d w x^=9 



2k 2 



4 V ; 48 v ; 



1 1 

2^2 



B A dC {3) A dC {3) (4) 



where is the dilaton field, g^ v is the metric, H = dB is the NSNS three- 
form field strength, F^> = dC^ 1 ' is the RR two-form field strength and = 
dC™ +H ACW is the RR four-form field strength. The system of iV DO-branes 
couples to all of these fields. 

The effective action of iV DO-branes for weak and slowly varying fields is 
the non-abelian U(N) Yang-Mills action plus the Chern-Simons action (for 
the bosonic part). For weak fields the action is gotten by dimensionally re- 
ducing the action of 9+1 dimensional U(N) Super Yang-Mills theory to 0+1 
dimensions 0. Up to a constant term it is 

S = -T (2tt/ s 2 ) 2 J dt Tr (\f^A , (5) 

where is the non-abelian U(N) field strength in the adjoint representation 
and To = (gX) -1 is the DO-brane mass. To write this action in terms of 
coordinate matrices X 1 , one has to use the dictionary 

Ai = -^X\ F oi = —-X i , Ftj = — - \ i\X\Xi] (6) 
2vr/ s 2 ' 2nl 2 s ' 3 (2tt/ 2 ) 2 L J K ' 

with i,j = 1,2, ...,9, giving 

S = T jdtTT (ipX* + l-L^^Xi^X^ (7) 
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To derive this it is necessary to gauge the A potential away, which is possible 
for a non-compact time. The Gauss constraint 



[X i ,X i ] = (8) 

persists, and the above action should be taken together with it ||. 

The assumed approximation that the fields F 0i and Fy in should be 
weak and slowly-varying corresponds for the coordinate matrices X 1 to 

x j <i, [x\x j ]<z:i 2 s , J*r<c\ [X\Xi\<£l s . (9) 

The condition X ! < 1 means that the DO-branes are in the non-relativistic 
limit. 

Recently it has become clear |2|, [| that even though the DO-brane world- 
volume is only one- dimensional, a multiple DO-brane system can also couple 
to the branes of higher dimension via the Chern-Simons action. The Chern- 
Simons action not only tells us how N DO-branes move in the weak background 
fields of Type IIA supergravity, but also what higher fields the DO-branes 
produce. Using this, one can see that it is possible to build the higher p-branes 
of Type IIA string theory out of DO-branes |TU], ^ |3| . The Chern-Simons action 
derived in [0, |, || for the coupling of N DO-branes to bulk RR C (1) and C (3) 
fields is 

Scs = T J dfTt^o + a^ + ^fc-w^ (10) 

The idea that a lower-dimensional object under the influence of higher- 
form RR fields may nucleate into spherically or cylindrically wrapped D-brane 
was proposed by Emparan fl2| for the case of fundamental string. Before that 



Callan and Maldacena [13|] constructed a D- or F-string as a BI soliton solution 
on the D-brane where the attached string appears as a spike on the brane. 

We first clarify the notion of the coupling of the D-brane world- volume RR 
current to the external form-field, in particular for wrapped D-branes. This 
will give us a better understanding of the similarity between the D2-brane 
current and the corresponding expression for a system of DO-branes. 

In general the D2-brane couples to the bulk RR field through the well- 
known CS coupling 

Scs = I C< 3 > = / C avn J^ p d?a (U) 



>CS ~ J ~ I >"'<•"' 

where J is a three form RR current 

r vp = e a ^ d a X»d p X v d y X p , (12) 

or in form notation 

*J^ P = dX» AdX u AdX p . (13) 
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One can introduce a charge corresponding to this current, such that it is a 
world- volume 2-form even though it has three space-time indices, same as the 
current 

= x^dpX v d 1 Xd , so that 

ZZ = tdwQZ ■ (14) 



or in form notation 



Q^p = (dX [lA A dX u ) X p] , such that 
where the exterior derivative is taken with respect to world-volume indices. 



This charge first appeared in the theory of bosonic relativistic membrane [14 



There it was interpreted as a topological charge of the membrane. The con- 
nection between RR charges and D-branes was discovered by Polchinski [15 . 

Instead of (|12D one can represent the current as a Poisson bracket with 
respect to the spatial world-volume coordinates. For a static membrane, the 
non-zero components of J pup are 

J 0ij = {X\X j } , (16) 

For a moving membrane the completely spatial components also appear. A 
convenient generalization is, in the static gauge 

J ijk = X i {X j ,X k ] . (17) 

The above discussion is for the ordinary D2-brane, and the coupling of the 
DO-brane matrix-mechanical system to the field is completely analogous 

Scs= ^kl dt Tr { c °^ xi > + °^ x ^ x ^ xk ) = S c ■ j dt (i8) 

The trace of the first term is zero identically, reflecting the fact that the total 
bare RR charge of the object that we are considering is zero, due to cancellation 
of the pieces with opposite orientation on the 2-sphere. One needs to expand 
the expression in powers of X, to effectively obtain the multipole expansion 
of C^ 3 \ Then, it is possible to integrate by parts, to get instead the coupling 
of the field strength to the charge || 

Scs= ^p I dt F ° i]k Tr xj ] xk= S F ' Q dt (19) 



In Section we will use this expression for the D2-brane dipole moment Q 
to compute dipole radiation. This is convenient because J can be easily 
obtained by time differentiation. 
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2.2 The spherical membrane 



In this section we briefly review the spherical D2-D0 brane configuration of 
type IIA string theory since our new solution of the system of N DO-branes 
presented in Section |3] uses the essential elements of this construction. The 
solution is equivalent to the spherical membrane solution of M(atrix) theory 

SI- 

We aim to construct a membrane with an S 2 geometry. We embed the 
S 2 in a three dimensional space spanned by the 123 directions. We take the 
ansatz 



Xi(t) 



-.Tinit), i = l,2,3 



(20) 



VW - 1 

where the N x N matrices T 1; T 2 , T 3 are the generators of the N dimensional 
irreducible representation of SU(2), with algebra 



[T. T,] 



and with the quadratic Casimir 

3 



N 2 



so that Tr(T- 



iV(iV 2 - 1) 



12 



(21) 



(22) 



For vanishing background fields the Hamiltonian isQ 

1 3 2 

1 -2 , a ( 2 2 , 2 2, 

2 l^ T i + Y V 1 ^ 



H 



NTn 



2 2 
r 2 r 3 



i=l 



where we have introduced the convenient parameter a 
the equations of motion 



(23) 

. This gives 



2/2 i 2\ 

ri = —a (r 2 + r 3 )ri 

2/2 2 \ 

r 2 = -a (r 1 + r 3 )r 2 
h = -^ 2 {r\ + r%)r 3 



(24) 



Let us estimate the physical size of the resulting object. For simplicity we take 
all radii to be equal to each other: r± = r 2 = r 3 = r. With this we have from 
(f2~0|) and (|22|) the physical radius of the membrane 



R 1 



xl + xl + xl 



I r z 



(25) 



where I is the N x N identity matrix. The formula (|25|) shows that the iV 
DO-branes are constrained to lie on an S 2 sphere of radius r. This condition 
is what determines the normalization of the ansatz (|20|) . 

Considering the Chern-Simons action ([19|) we clearly see that the coupling 
of this system to -F0123 is non-vanishing. This means that the spherical mem- 
brane solution has a D2-brane dipole moment, and it is thus appropriate to 

3 This system is otherwise known as + 1 dimensional classical SU(2) YM mechanics 0]. 



7 



recognize the spherical membrane solution as a bound state of a spherical 
D2-brane and N DO-branes M. 

The equations of motion ( 24] ) reduce to the equation r = -2a 2 r 3 . Thus, it 



is clearly not possible to have a static solution^]. Instead, the solution is [16 



r(t) = R sn(R Q at + 0) (26) 

where sn is the sinus amplitudinis of Jacobi defined by 

ry i 

y = sn(x) -vv- x = I —f===dz (27) 



z 



The solution ( |26| ) is thus an S 2 sphere with the radius r(t) oscillating between 
Ro and —Rq with a period of ^J= Clearly the spherical membrane will 
be classically point-like at the nodes of the sinus. Thus the membrane solution 
will break down after a finite amount of time, since the classical solution may 
not be valid at substringy distances, and possibly decay into a Schwarzschild 
black hole 0. 

The conditions @ for the spherical membrane solution translate into 



\r(t)\ < l s VN, \r(t)\ < 1, \f(t)\ < C (28) 

where we have used the fact that the matrices T are of order N. Since we also 
require \r(t)\ ^> l s we must have N 3> 1. Thus it is necessary to have a large 
amount of DO-branes to build a macroscopic spherical membrane. In fact, the 
maximal allowed size with the use of N DO-branes is proportional to y/N. 



3 The rotating ellipsoidal membrane 

The spherical membrane solution reviewed in Section |2]2| is not a stable object 
of string theory. It is classically unstable under small perturbations, and also 
reaches a point-like singularity in finite time. This makes it possible for it 
to collapse into a near-extremal Schwarzschild black hole. In this section we 
present instead another kind of solution to the system of iV DO-branes which 
is free of singularities. We also show in Section [| that it is a stable object, in 
the sense of stability under small perturbations of initial conditions. Moreover 
it has interesting physical properties due to its dynamical nature, like the 
radiation of various SUGRA fields (see Section |5|). The basic idea in the 
construction is that the attractive force of tension should be cancelled by the 
centrifugal repulsion force. The motion is at all times transverse, and cannot 
be gauged away by coordinate reparametrization invariance on the membrane. 

We now construct a rotating ellipsoidal membrane, viewed as a collection 
of DO-branes, in such a manner as to cancel the attractive force of tension with 

4 Myers B considers the system of N DO-branes in a constant external 4-form RR field 
strength, as in (|l9|). Then indeed there is a static solution, where the DO-branes polarize 
and arrange into a static spherical configuration. 
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the centrifugal repulsion force. For that we need to take the basic configura- 
tion ( pop of the non-commutative fuzzy sphere in the 135 directions, and set 
it to rotate in the transverse space along three different axis, i.e. in the 12, 
34 and 56 planes. We thus use a total of 6 space dimensions to embed our 
D-brane system. The corresponding ansatz is^ 

= V^T Tiri(t) ' XM = 7W=T TirM • 

X ^ = 7rhf* T3{t) • Xi(t) = 7¥hf 2T ' {t) ■ 

Xs(t) = TW^j" 1 ' 3 ^ ' Mt) = 7W^r T3r ° l - t) ' (29) 



In this new ansatz we are still using the SU(2) matrix structure from ( |20D 
such that the coordinate matrices are proportional to the SU(2) generators in 
pairs. We interpret this as a rotation because one could make a rotation in 
e.g. the 12 plane which makes one of the components vanish, say X 2 , while 
the other one gets a radius r[ = \Jr\ + r\. This is possible exactly because 
both are proportional to the same matrix TV The end result is that at any 
point in time one can choose a coordinate system in which the object spans 
only three space dimensions. 

Substituting the ansatz into (0) gives the Hamiltonian 




( r 2 + r 2)( r 2 + r 2) 



+ (r? + r 2 2 )(r 5 2 + r 6 2 ) + (r 3 2 + r 4 2 )(r 5 2 + r 6 2 )j J (30) 
The corresponding equations of motion are 

2i 2 1 2 1 2 1 2\ •■ 2/ 2 , 2 , 2 , 2\ 

ri = -a (r 3 + r 4 + r 5 + r e ) r x , r 2 = -a (r 3 + r 4 + r 5 + r 6 J r 2 , 

2/ 2 , 2 , 2 i 2\ •■ 2/ 2 , 2 , 2 , 2\ 

r 5 = -a 2 (r\ + r 2 + r 2 + r 2 ) r 5 , f 6 = -a 2 (r\ + r 2 + r 2 + r 2 ) r 6 . (31) 

We have found the special solution to these equations describing a rotating 
ellipsoidal membrane with three distinct principle radii Ri, i? 2 and R3 

ri(t) = Rx cos(co>ii + 0x) , r 2 (t) = -Ri sin(tuit + 0i) , 
r 3 (t) = R 2 cos(u 2 t + 02 ) , r 4 (£) = /2 2 sin(cu 2 )f: + 2 ) , 
r 5 (t) = i? 3 cos^t + 3 ) , r 6 (t) = i? 3 sin(u;3t + 3 ) . (32) 

This particular functional form of the solution ensures that the highly non- 
linear equations for any of the components Ti are reduced to a harmonic oscil- 
lator. The solution fl3"2] ) keeps r\ + r\ = R\ , r\ + r\ = R\ and r\ + r\ = R\ 

5 A less general ansatz was considered in |l7j] to model a spherical rotating membrane, 
see also pq|. Moreover, a generic rotation-ansatz was proposed in p9|. 
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fixed which allows us to say that the object described by ( p2[) rotates in six 
spatial dimensions as a whole without changing its basic shape. 

Using the equations of motion (|31~D, the three angular velocities are deter- 
mined by the radii, and do not necessarily have to coincide: 



out = a^/Rl + R 2 3 , u 2 = ayJRl + R 2 3 , u 3 = a^/Rj + R\ . (33) 

This dependence of the angular frequency on the radii we interpret as if the 
repulsive force of rotation has to be balanced with the attractive force of tension 
in order for ( |32"D to be a solution. Thus the radii R\, R 2 and R% parameterize 
(p^) along with the three phases fa, to produce altogether a six parameter 
family of solutions. Unless explicitly stated otherwise, we set fa = in what 
follows. 

In order to exhibit the properties of the solution (|32|) we compute the 
components of the angular momentum 



Mij = Tr 



(34) 



where IP = T X\ As expected, the only non-zero components are 



M l2 = ^NT Q u 1 R\ , M 34 = ^NT lu 2 R 2 2 , M 56 = ^NT co 3 R 2 . (35) 

The angular momenta M 12 , M 34 and M 56 correspond to rotations in the 12, 34 
and 56 planes respectively. Their values (|35"D fit with the interpretation of the 
solution being N DO-branes rotating as an ellipsoidal membrane in that they 
are time-independent due to conservation law and proportional to NTquoiR 2 . 
The conditions (§) on the solution (|32]) give 

UiRi < 1, RiRj < Nl 2 , uj 2 Ri < l; 1 , UiRiRj < Nl s , i ^ j (36) 

Using (^) this is seen to be equivalent to 

RiRj <C Nl 2 , RiR 2 <^N 2 ll, i^j (37) 

Since we also require the membrane to have a size larger than the string length 
l s , we must have N ^> 1, just as for the spherical membrane reviewed in 
Section |2T2"|. In the special case Ri = R 2 = R 3 = R the conditions ([J7]) reduce 



to R < VNl 



4 Stability analysis 

In this section we consider the problem of stability of our solution (^) with 
respect to an arbitrary change of initial conditions. In principle, there may 
be other possible sources of instability, classical or quantum. The quantum, 
or rather quasi-classical, loss of energy associated with radiation of various 
supergravity fields is computed in Section |]. 
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In Hamiltonian systems there is no naturally defined positive-definite metric 
on the phase space, therefore so-called Lyapunov exponents are only meaning- 
ful when measured on time-slices of a periodic trajectory, as used by Poincare. 
Thus in what follows we should take ^1,2,3 to be multiples of a lowest common 
frequency u , such that the trajectory is periodic with T = Our formulas 
are applicable to just such a situation, but when it comes to measuring the 
exponents we take the simplest case ojq = uj\ = u 2 = uj 3 , in which the ellipsoid 
becomes a sphere. 

Let us consider linear perturbations to the equations of motion fl3"T|), which 
are essential in all kinds of stability analysis. These explicitly time-dependent, 
or else called non-autonomous, equations should be understood to contain the 
original solutions fl32|) as time-dependent external fields: 



a~ 2 5r 1 = —6ri(Rl + - 2r 1 (r 3 5r 3 + r 4 5r 4 + r 5 <5r 5 + r 6 <5r 6 ) 
a~ 2 5r 2 = —5r 2 {R\ + R\) - 2r 2 (r 3 5r 3 + r 4 5r 4 + r 5 5r 5 + r 6 5r e ) 
a~ 2 8r 3 = -8r 3 {R\ + R 2 3 ) - 2r 3 (r 1 8r 1 + r 2 5r 2 + r 5 5r 5 + r 6 5r 6 ) 
a~ 2 5r 4 = -5r A (R\ + R\) - 2r 4 (r 1 5r 1 + r 2 5r 2 + r 5 <5r 5 + r 6 <5r 6 ) 
a~ 2 5r 5 = —8r 5 (Rl + R%) - 2r 5 (nSri + r 2 5r 2 + r 3 5r 3 + r 4 5r 4 ) 
a~ 2 5r 6 = -5r & (R\ + R\) - 2r 6 (ri5r a + r 2 5r 2 + r 3 5r 3 + r 4 5r 4 ) (38) 

These describe the time evolution of 5r(t), the linear deviation from a given so- 
lution r(t). We could rewrite these equations in first-order form, and introduce 
the 12-component perturbation vector 

a(t) = (5n(t), 5h(t), 6r 2 (t), 8r 2 {t), 5r 6 (t), 5r 6 (t)) (39) 

We have then the first-order linear system of differential equations 

a(t) = H(t)a(t) (40) 

where H(t) is a 12 x 12 real matrix, put together from a 6 x 6 matrix to be read 
off RHS of ( |3"8"D and a 6 x 6 unit matrix in the appropriate order of components. 
This matrix is periodic with the period equal to that of the original solution. 
Time ordered exponentiation gives us the formal expression for the complete 
evolution S-matrix 

S(t)= : feapjf H(r)dr) : , (41) 

so that S(t) = H(t). We can now find the solution vector by applying the 
S-matrix to the initial vector: 

a(t) = S(t)a(0) . (42) 

In order to investigate the stability of the original periodic trajectory we need 
to find out whether small perturbations grow. This information is encoded in 
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a single marix S(T ) which gives the evolution of the perturbation vector over 
one period. This matrix is symplectic, and has unit determinant in order to 
preserve the Liouville measure on the phase space. 

The original solution is parameterized by Ri,R 2 , R3 and 4>i, fa, fa. A small 
variation of either of these parameters produces a new solution. Together, they 
give 6 different zero modes of the system. Technically, zero modes are produced 
by formally differentiating the solutions (|32| ) with respect to a parameter. 
Physically, this is very familiar phenomenon, e.g. the free scale parameter in 
the instanton of YM gives a zero mode in exactly the same fashion as here. 
Zero modes correspond to unit eigenvalues of the S-matrix. The zero modes 
for d/dfa, i = 1, 2, 3, are periodic since for example 

d 

5r iW = ^T r iW = -Ri sinM + fa) (43) 
ofa 

The initial conditions for d/d<f)i, i = 1, 2, 3 are 

a 01 (O) = (0,-ujxRuRx, 0,0, 0,0, 0,0, 0,0,0) 
a 02 (O) = (0,0,0,0,0,-^,^,0,0,0,0,0) 
a 03 (O) = (0,0,0,0,0,0,0,0,0,-^3,^3,0) (44) 

On the other hand, the d/dRi initial vectors are 

a fll = l,0,0,^,0,0,0,« 2 ^^,0,0,0,a 2 ^^ 
a/? 2 (0) = (0,0,0,a 2 ^i,l,0,0,^ 2 ,0,0,0,a 2i? " /,>:; 



a fl3 (0) = ( ,0,0,a 2 ^,0,0,0,a 2 ^, 1,0,0,^ • (45) 



These three modes correspond to nearby trajectories of the same form as 
and have a slightly different period. Consequently the trajectory remains in 
the vicinity, but gets either ahead or behind of the original one linearly in time 
due to the difference in period. The direction of this linearly growing deviation 
is, not surprisingly a combination of a^.. For example after k periods 

a Rl (A>T ) = S(fc-r )a(0) = S(T ) fe a(0) = ai?1 (0) + A;-T [a^(0) + a^(0)] (46) 

Now that we have singled out the zero modes, we can proceed to find the 
other 6 eigenvalues by integrating the equations numerically and discarding 
the unit eigenvalues. The remaining 6 eigenvalues were found, together with 
their eigenvectors, and are exhibited in Table 1. 

It is a general fact that the characteristic polynomial of a symplectic matrix 
is reflexive, and therefore the eigenvalues come in one of the three possibilities 
12T 



a) Ai 5 2 are Real Ai • A2 = 1 
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\ 

A 


arg A 


-.11253 - .99364 % 
— .LLzoo + .yyoo4 i 


-96.461° 

A fil O 

yo.4ol 


-.11253 - .99364 % 
-.11253 + .99364 % 


-96.461° 
96.461° 


-.95006 - .31206 i 
-.95006 + .31206 i 


-161.81° 
161.81° 



Table 1: Eigenvalues. To the available accuracy O(10~ 5 ) all of the eigenvalues 
are of unit length. In the second column the value of argA is given in degrees, 
which has the meaning of the angle of rotation between the corresponding two 
eigenvectors per period T . 



13) Ai, 2 are Unitary Ai = e id A 2 = e~ ie X 1 ■ X 2 = 1 

7) when two pairs of the second kind collide on the unit circle, they can go 
off and form a quartic arrangement: 



Ai 

A 3 



A 2 
A 4 



Ae- 
1 



Ai ■ A2 • A3 • A4 — 1 . 




Figure 1: The possible location of eigenvalues. 



The system is stable if and only if it has pairs of eigenvalues of the second 
kind (P) only. Eigenvalues outside of the unit circle (a), (7) correspond to 
eigenvectors that grow exponentially with time (Figure 1). Conversely, eigen- 
values on the unit circle correspond to a pair of real vectors which are rotated 
into each other without deformation. The angle of rotation per period Tq is 
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equal to argA and quoted in the table as well. The coinciding pairs of eigen- 
values in the table are accidental due to symmetry, and we dont expect it to 
happen for general Ri. 

The conclusion is that using a combination of analytical and numerical 
methods we were able to demonstrate stability of the specific trajectories of 
the equations of motion ([31]). We stress that in general, classical YM mechanics 
exhibits strong instability of periodic trajectories. In particular, the pure SU(2) 
system, eqns fl24D was conjectured || to have no stable periodic trajectories, 
and may even be an ergodic, k-mixing system. 



5 Radiation 

5.1 Ramond-Ramond one-form radiation 

In this section we compute the C^-field radiation from the rotating ellip- 
soidal membrane. DO-branes are electrically charged under this field, but as 
we will see, the dipole moment of the system vanishes and we have to look at 
quadrupole order to find the radiated wave. 

The dynamics of the C^-field is given by the action 



l - J d 10 x^(F^) 2 + T j dtTr(c + Cyr) (47) 



2 

Choosing the Lorenz gauge 

d^C = (48) 
we have from ([TJ) the wave equation 

d x d x C t = 2k 2 T q Tr{X l )5{x) (49) 

In the region far from the source we have the plane wave solution 

C,j,(x,t) = C^exp(ik\x x ) + c.c. (50) 

From flU) and flU) we get 

fc M fcf = 0, = (51) 

The energy-momentum tensor is 

2n% u = l -F x F uX - ± %U F 2 (52) 

Averaging over a region much larger than uo^ 1 we get the energy-momentum 
tensor of the plane wave 

(T^u) = i^k^ u C x C* x (53) 
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We set k° = uj and k = ux/r in what follows in order to make the local approx- 
imation of the spherical wave as a plane wave. Thus, we are only considering 
the radiation contribution at the particular angular frequency uj. 

The source for Cj in ( |49"D is clearly zero, so there is not any dipole radiation 
for the rotating ellipsoidal membrane. Instead, we must analyze the differences 
in propagation times due to the fact that the object is extended in space, since 
this gives the quadrupole radiation. In order to derive the quadrupole source, 
we use electrodynamical analogy and consider instead N particles of equal 
charge at positions x a , a = 1...N. For these, the wave equation is 

N 

d\d x Ci = 2k 2 T ~ Za)e- Wt (54) 

a=l 

where the angular frequency uj is the considered oscillation frequency. Using 
the propagator ( |104j ), we get 

o 2rn 3 r N 
V ; 105fi 8 r 4 J j^i 

9 K 2rn 3 N 

Making the identification Y2 a =i x a x a ~~ y 

Tr(J<W) we get 

Ci{s ' t] = ~*S^ e ^ e ^ ( Tr( ^ } ~ lk] Tl{ ^ XJ) ) + cx - (56) 

In general Ti{X' l X j ) has both a symmetric part Q ij = i(Tr(X i X J ')+Tr(X- J 'X i )) 
and an antisymmetric part |(Tr(X*X J ) — Tr(X^X % )). The symmetric part 
corresponds to electric quadrupole radiation and the antisymmetric part to 
magnetic dipole radiation. The antisymmetric part is proportional to the 
angular momentum tensor M %: > which is a constant of motion for our particular 
solution. Thus, for our solution the antisymmetric part does not contribute to 
the radiation, and we shall therefore set the antisymmetric part to zero in the 
following. Using ( pTj ) and setting Tr(X J ) = we finally have 

e > = - (57) 

where k = k/u. From ( |51"D we have that ujCq + k l Ci = 0, from which we get 

C Q = -k % Ci = t^-^^rkhQ^ (58) 
105^8 r 4 ^ K ' 

Therefore, from (0), Q58D and we get that the power radiated per solid 
angle is 

dfl* N 1 2k 2 a 



2k 2 T ( 



105 2 ^ 



° uj 10 



5iikjk m kikjkik m Q ^Q (59) 



15 



Integrating (|59|) over the S 8 sphere, we get the total radiated power 



1 ?K 2 T 2 



3 2 5 2 7 2 n n» 



q ij q: 3 -\\q i 



(60) 



Note the similarity between the expressions ( p9| ) and (^) and the correspond- 
ing expressions (|108|) and ( [111| ) for the gravitational quadrupole radiation. 



We now need to compute the quadrupole moment Q 1 ^ for the rotating 
ellipsoidal membrane. One can check that it splits into a block-diagonal form 
with three two-by-two matrices corresponding to the three coordinate pairs 12, 
34 and 56, with all other entries in being zero. The 12 block is given by 

Qn Q 12 \ = 1 2 f -sin 2^* cos2^ \ 
Q21 Q22 J 6 1 V cos 2uJ ^ sm 2w i* / 

The 34 and 56 blocks are similarly given in terms of the two other frequencies 
and radii. One can see from ([BID that the quadrupole radiation has double the 
frequency of the rotation frequency. Thus, the rotating ellipsoidal membrane 
emits quadrupole radiation at frequencies 2u 1 , 2u 2 and 2co 3 . 

Considering the fourier component at frequency 2u±, we get that the non- 
zero components of Q 1 ^ are 



From formula ((6^) we then get 



Pl = 1 2K " lo u\ 2 R* (63) 
3 4 5 2 7 2 11 Q 8 

The total power emitted via quadrupole radiation in the C^-field is therefore 

p = 3*£m 2K T 8 T ° 2 H2R{ + ul2Rt + ufRt) (64) 

As one can see from flM] ) the total power has the dependence uj 12 R on the fre- 
quencies and radii. Instead, electric quadrupole radiation in 3+1 dimensions 
has the dependence uo 6 R 4 . But, these two cases are in perfect correspondence 
with each other since the propagator (|104j) in 9+1 dimensions has an a; 3 depen- 



dence in the leading 1/r 4 term. This gives an extra u 6 for all types of radiation 
when comparing with 3+1 dimensions. For example, scalar radiation has the 
dependence u 2 on the frequency in 3+1 dimensions and u; 8 in 9+1 dimensions. 
That radiation in 9+1 is suppressed by the factor u; 6 in comparison with 3+1 
dimensions, can be understood as the fact that there are 6 more directions 
available for a particle to propagate in. 
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5.2 Ramond-Ramond three- form radiation 



We would like to compute the radiation due to the coupling of the membrane 
to the RR potential. At this point we make use of the couplings we derived 
in Section [2.1|. The dynamics of the -field is governed by the action 



Sew 



2k 2 



d 10 xV-g 



2-4! 



Xpvp 



pXpvp 



+T 



2nli 



Irk 



dtTr(C 0lJ [X\X3}+C l3k [X\Xi]X 



where 



?< 4 > 

Xpup 



(3) 



[A pup] 



(65) 
(66) 

(67) 



where we now have to impose anti-symmetrization with respect to ijk indices. 

We should be working in Lorenz gauge, in order to ensure the explicit 
transversality of the radiated wave: 



We denote the current by J ljfc with a convenient normalization 

2nli V / 



pup 







This gives the wave equation 



d x d x C ljk = 2K 2 T J^ k 5(x) 
In the region far from the source, we have the plane wave solution 

Cp; V p{x,t) = C^ p exp(ik x x x ) + C* up exp(-ik x x x ) 
with the conditions 



kpk^ — 0, k^C^p 







(68) 



(69) 



(70) 



(71) 



In order to compute the energy carried away by this wave we will need the 
energy-momentum tensor: 



2k 2 T, 



pu 



1 1 

2 3! 



p x P<?p 



uXpa 



2^2-4!" 



(72) 



We now set k° = u and k = ujx/t. We see from ([H]) that fields that contain 
time-components can be found by using the gauge condition C jk = —kiC^. 
The field strength is gotten by appropriate differentiation of the plane wave 
solution flTDD, because in the radiation zone the leading order derivative always 
comes from the exponential: 



Fijkl 

Fojki 



4dfiC 



[i^jkl] = 



k{Cjki 



kjCiki kkCm k]C 



H^jki 



e Mr~t) + c _ c _ 



h-,, — % h,f) '., — h h,H,- —h %H,, 
^jkl rv m rv l^ J ]km rh rn nj k K - / ljm n "in rh j K - J klm 



iuj(i — i) 



+ c.c. (73) 
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After averaging the energy-momentum tensor over a space-time region of size 

uj- 1 

2/« 2 (7]iy) = k fl k 1/ —C fll/p C^ up = kfjzv— ( ' , t i-( 'jjj. — 3k l Cijkk m C m j k (74) 

Note that this form of (T M1/ ) was to be expected since the only Lorentz invariant 
scalar one can construct from C pup is C pup C* up . The F 2 term in (|72|) vanishes 
just like in ordinary Maxwell plane wave, as can be ascertained by direct 
contraction from ([731). 

The propagator in 9+1 dimensions, including the normalization was derived 
in the Appendix [| and gives the following form for the C-field in the radiation 
zone: 

C ijk (f,t) = 2K 2 T ^-^e^ r -^ + c.c. , (75) 
3 V ' ' ° 105fi 8 r 4 V ' 

from which we read off the polarization 

We are now ready to compute the energy flux in the out direction 



dP 2 t 2k 2 T 8 1 
dOT " 2K Toih ~ 105*nf 3! 



J i]k J* jk -3hk m J lij JL 



(77) 



'8 iuJ "8 

For our particular system the non-zero components of Q are 

Q135 = -NaRiR 2 R 3 cos u)\t cos u 2 t cos oo 3 t 
3 

Q246 — -N01R1R2R1, sincjit sinco> 2 £ sinc^t (78) 

together with Q im , Q im , Q235 , Q236 , Q145 , ^245- We see that it is possible 
to produce higher harmonics by multiplying the above trigonometries. For 
generic values of u>i one has four different frequencies: 

Co>l + Co> 2 + Co>3 , + ^2 — ^3 , ~ ^2 + ^3, and — U)\ + U 2 + UJ3 ■ 

The computations are simplest for the highest mode, and in any case this is 
the one that carries the most power. We pick out this mode u = io\ + uj 2 + u 3 , 
and also set R 3 = R\R 2 R?j 



Q135 = ^NaR 3 e^ , Q 235 = Qi 45 = Qise = ^NalPe™ 

Q 2 46 = —Natfe™ , Q 245 = Que = Q236 = — 2 \NaR 3 e' kjt (79) 

From these it is easy to obtain the corresponding J's by differentiation with 
respect to time. Substituting the fields into ([77D we finally get the power 
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radiated into unit solid angle as a function of direction, here parametrized by 
direction vector k 



dP 



?k 2 T 2 1 



105 2 Qi 



1 6 



This dependence on the direction is typical of rotating sources. The polariza- 
tion of the wave (|79|) is such that when looking at the wave in a completely 
transverse direction (the 789 hyperplane), the wave is "circularly" polarized. In 
the 123456 hyperplane, the wave is "linearly" polarized and has half the power. 
In fact, we can rewrite the directional dependence in (|80|) in terms of the angle 
between k and the 789 hyperplane. The result is simply [1 — sin 2 #/2]. So 
indeed this is similar to the elliptically polarized field produced by a rotating, 
rather than oscillating dipole. 

The total power emitted at the frequency to = to\ + 0J2 + 1^3 is, after inte- 
grating over S 8 

P= ± ^J^^0 N 2 a 2 R 6 (gl) 



12 105 2 a 



5.3 Gravitational radiation 

In this section we compute the power radiated from the rotating ellipsoidal 
membrane in the form of gravitational field. In the Appendix |6| we derived the 
formula (|1 1 1|) giving the radiated power for a given energy-momentum tensor, 
in terms of its fourier components. For the rotating ellipsoidal membrane the 
energy momentum tensor is 

T M „ = T (2vrZ 2 ) 2 Tr (f^F* - j^T 2 ) (82) 

We quote the Too component first, since it is the only one that has an interesting 
time-independent piece, namely the total energy 

Too = ^ATTo 5> 2 F 2 (83) 

i=l 

The spatial non-zero components are split into three blocks: 12, 34 and 56. 
Each one is oscillating at its own frequency, and since we are interested in the 
fourier decomposition of T^, we can start with the only block that contributes 
to radiation at frequency 2co>i 

( Tn T12 \ _ NT p2 2 / cos2wit sin2^it \ , . 

V T 21 T 22 J' 6 11 V sm2^it -cos2c.it J [M) 

Note that the energy-momentum tensor naturally has twice the frequency of 
the underlying object. This is a general feature related to the fact that there 
is only one kind gravitational charge. 
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Other blocks are of the same form, so for non-coincident frequencies only 
one block corresponds to each frequency: the 34 block to frequency 2ui2 and 
56 block to 20J3. Note also that it is traceless, and so the power will be given 
by the first term alone in the formula fllllj) which we rewrite for convenience 
below: 

^ s^i-J 1 *^-^'') (85) 

Finally, the total power emitted at frequency 2oj\ is 

The exact same formula holds for radiation at the other frequencies. There can 
be no interference effect between different frequencies, thus the total emitted 
power is given by summing the conributions 

96 9 K 2i\r2rp2 

" \u{ 2 Rt + u?RZ + u?I® (87) 



3 2 -5-7-ll n 



We note that each term is of the form P ~ k 2 u; 8 (/u; 2 ) 2 where the moment of 
inertia I ~ NT Q R 2 . Overall, the radiation is extremely suppressed by twelve 
powers of the frequency. In Appendix || we explain that this differs from 3 + 1 
dimensional power of six, solely because of extra six powers coming from the 
propagator in 9 + 1 dimensions. 



6 Discussion 



We have shown that our rotating ellipsoidal membrane solution is classically 
stable. However the system constantly loses energy due to the quasi-classical 
radiation of various supergravity waves. This immediately poses the question 
of what will happen to the rotating ellipsoidal membrane after all available 
kinetic energy has been radiated. There are two possibilities: Either it will 
decay into iV free DO-branes or it will eventually become completely stable. 
We stress that the energy is positive, and so existence of a ground state is not 
obvious. If the system is to be quantum-mechanically stable, the size must be 
of order l s , since at larger length-scales the classical approximation is valid and 
it does radiate energy. Therefore, if indeed this scenario is actually realized, 
there would be a particle-like ground state with positive energy. 

In order to clarify what we mean by this fundamental particle ground state, 
let us employ the analogy with the archetypical example of quantum mechan- 
ics: the hydrogen atom. In this paper we found a rotating solution, similar to 
the motion of the electron in the central electrostatic field, however the elec- 
tron should lose energy and emit photons, eventually falling onto the nucleus. 
This is the paradox which was resolved by the Bohr-Sommerfeld quantization 
of the atom. In this paper we computed a similar instability due to emission 
of various supergravity waves. 
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In the final quantum theory, the electron indeed emits quanta of light, 
corresponding to transitions between discrete energy levels. But there exists a 
ground state which is completely stable. 

The Bohr-Sommerfeld quantization scheme tells us that the angular mo- 
mentum should be quantized in units of the Planck constant. We thus have 

M l2 = ^NT uo 1 R\ = L l , 

M 34 = ^NT u 2 R 2 2 = L 2 , 

M 56 = ^NT cu 3 R 2 3 = L 3 , (88) 

where Li are integer values of the projection of the angular momentum. Roughly 
speaking, for L = L\ = L 2 = L 3 , the allowed frequency spectrum is 

"i = 2T C-iiVT,(U)" (89) 

so that the leading behaviour is 

UL ~ ' (90) 

from this the quantized radius and energy are 

R L ~</(hLls and E ~ Luo ~ L^L^- (91) 

One question that we might ask is whether the angular momentum of the 
hypothetical ground state would be zero. The likely answer, by analogy with 
hydrogen atom, is positive. 
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A Gravitational radiation in 9+1 dimensions 

We consider the region far away from the source. We have the metric 

= r)ixv + V ( 92 ) 

where r)^ = diag(— 1, +1, +1, +1) and \h^ u \ <^ 1. Indices are always lowered 
and raised with r/^ and i]^. The Einstein equations are 

~ IjP^^f Rpa = T\iv + t^u (93) 
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where the left-hand side are first order in h^ v and is the energy-momentum 
tensor of the gravitational field. The linearized Einstein equations are 

~ \^ pa R { p} = T» v (94) 

Choosing the harmonic gauge 

& Ku = \d v h\ (95) 



We have from (|94f) the linearized Einstein equations 



1 " i.\ ;. _ o..2 ( rr 1 , A 



= 2k 2 IT^ - -v^T\j (96) 

Away from the source it is simply the wave equation 

d x d x V = (97) 
We consider a plane gravitational wave 

h^ u = exp(ik x x x ) + e*„ exp(-ik x x x ) (98) 
From EOM fl97|) and the gauge condition (|95|) we get 

fc^" = 0, fc'V = ^,e A A (99) 

We also have that e^ u = e Vil . Using the additional residual gauge freedom one 
can see that there are 35 independent components of e^ v . To second order in 
hfj, u the energy-momentum tensor of the gravitional field is 

*V = flg - \v,uV pa R?J (ioo) 

where we used that R$ = from fl9~7|). Averaging over a region much larger 
than we compute 



< i2g >= ( eV * e Ap - ^|e A A | 2 ) (101) 

Using this with fl9~9"|) and ( |100| ) we get the energy-momentum tensor of the 
gravitational wave as 

< V >= 2^M* ( eAp * e Ap - ^ \l 2 ) (102) 

In a 9+1 dimensional Minkowski space the solution to the wave equation with 
an oscillating 5-function source 

d x d x f{r,t) = 5{x)e- iu (103) 
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IS 

fir t) = -i-^-L^r ( 1+ ^_J5 15i_\ t , m 

n , } 105ft 8 r 4 \ l + ujr (ujr) 2 (urf ) C [ " } 

when demanding that only the outgoing wave is present in the solution^]. Using 
this, we can solve ( |9"ED for a single frequency uj as 

Ku(x, t) = -^-^e^e'^ J dV (t^x 1 , uj) - ~^T A A (f , u )j e***l r 

+c.c. (105) 
for r large. Setting k° = uj and k = ujx/r we have from ( |9~8"D that 

with 

T^fto;) = y d 9 f T^(x\u)e- lt *' (107) 
Thus, from ( |102| ) and ( |106| ) we get that the power emitted per solid angle is 

(1P r% < t ai >= -LwV (e^e Xp - ±|e A ' 2 



8k 2 / 1 

cu 8 A^ m (fc)T^T /m (105 



105 2 ftf~ \f V 
8k 2 



105 2 ^ 
where k = k/uj and 

7.... l 1„„ l. . 

A-ijim(k) = Sn5j m 2kjk m 5n-\- —kikjkik m —Sij5i rn -\-—kik m Sij-\-—kikj5i m (109) 

o o o o 

After averaging the above over the 8-sphere 
If - 7 

— / dfls^ijim{k) = - - — [895i/5 im - lOSijdim + 6 im Sji] (110) 
we finally get that the power emitted at the single frequency uj is 

= 2 ^ ( T n* T _ hr.A (in) 

Here we would like to remark on the overall dependence of the radiated power 
on frequency uj. The spatial components of the energy-momentum tensor above 
( |102j ) are typically proportional to Iuj 2 , and so power scales as Pglal ~ uj 12 - 
This is in contrast to 3 + 1 dimensions where the corresponding quantity scales 
as Pgrlv ~ uj 6 - However this is more of an artifact of the propagator ( 104|) as 
even I = scalar wave is radiated in 10 dimensions as Pglaiar ~ 008 > whereas in 
3 + ldim it's P^li ar ~ uj 2 ■ Thus the extra powers of uj in ( |111 ) are completely 
due to the a; 3 behavior in the propagator (|104| ). 

6 The solution in a d + 1 dimensional space-time is f(r,t) oc r~( d ~ 2 ^ 2 H^ > _ 2 y 2 (ujr)e' l ' Ajt , 



where Hn is the Henkel function of the first kind of order n. 
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